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Electronic structure calculations using quantum Monte Carlo (QMC) methods conclusively show 
that pure CaBg is a narrow-gap semiconductor with an X-point gap ~ 1.3 eV. This should put to 
rest controversies as to whether the compound is a semimetal or a semiconductor, as also theories 
suggesting it to be an excitonic insulator. Our extensive Hartree-Fock (HF) and density functional 
theory (DFT) calculations support the view that the hexaboride ferromagnetism can be induced by 
contamination with magnetic element atoms like Fe and that presence of La is not essential for this. 
La impurity, however, gives rise to a metallic state with a small electron-like Fermi surface as seen 
in ARPES experiments. 

PACS numbers: PACS: 71.20.-b, 71.15.Qe, 02.70.Ss 
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Recently, a great deal of interest has been generated 
in divalent hexaboride CaBg mainly due to observation 
of unusual ferromagnetism (FM) in the system suppos- 
edly with Lanthanum doping Q . The unusual feature of 
this ferromagnetic state is that it has very small moment 
(~ 0.07 /is per La atom) and yet has very high Curie 
temperature, T c ~ 600K. Theoretical models with un- 
conventional many-body mechanisms were proposed as 
explanation for this FM in absence of d- or /- electrons 
in the system. For example, it was suggested 0, that 
doped hexaboride forms the magnetically ordered ground 
state of a dilute electron gas while Zhitomirsky et al. and 
other workers 3| suggested it to be an excitonic insu- 
lator. Central to the idea of excitonic insulator is the 
assumption of a small band overlap of top of the B p va- 
lence states and bottom of the Ca d conduction states. 
This band overlap was suggested by early density func- 
tional theory (DFT) band structure calculations It is 
well-known by now that DFT in local density approxima- 
tion seriously underestimates the band gap, hence models 
based on a small band overlap in these calculations are 
somewhat questionable. Going beyond DFT, the pertur- 
bational GW method band structures, which incorporate 
self-energy corrections, have been calculated for the sys- 
tem. Though one of them predicts a large band gap at 
the X-point of the Brillouin zone (BZ) of this simple cu- 
bic system , another shows a small overlap . A very 
recent calculation using weighted density approximation 
predicts a band gap of ~ 0.8 eV at the X-point These 
conflicting results clearly point out that our understand- 
ing of basic properties of this interesting material is still 
rather limited. 



siderable controversy regarding metallic/semiconducting 
character of CaBg in the experimental results, too. Semi- 
metallic nature is supported by observation of metallic 
resistivity |8j, |9| , the de Haas- van Alphen effect H03, 
and Shubnikov-de Haas results ^]|. On the other hand, 
NMR 0, angle resolved photoemission spectroscopy 
(ARPES) an d thermopower experiments SU S~ 
gest that CaB 6 is a narrow-gap semiconductor. Some 
recent high resolution ARPES |l5j and combined trans- 
port, optical and tunneling measurements [l6| very 
strongly suggest that pure CaB 6 is a semiconductor, 
while metallic character and FM in the system originate 
from boron-related defects. 

While the pure system is poorly understood, there is 
very little understanding of the origin of FM in the sys- 
tem. Authors in ref. [ljj proposed a scenario for FM 
in La-dopcd samples over a small concentration range 
by studying formation energies and associated moments 
of various defects. However, given the large formation 
energy of Bg octahedron vacancy they found, and obser- 
vation of FM in systems even without La-doping, make 
their proposal less likely. Various experimental groups 
have suggested, and there is growing consensus on this 
point, that FM in CaB@ is an extrinsic effect related to 
magnetic impurities such as Fe and Ni 0, 0, , prob- 
ably coming from B used in sample preparation. The 
most convincing argument for this is provided by the ex- 
periments of Cho et. al |16| . 

Hoping to clarify some of these issues, we did a de- 
tailed first principles electronic structure calculation of 
pure CaBg using Hartree-Fock (HF), DFT and quantum 
Monte Carlo (QMC) methods, and also did an analysis 
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fects include substitutional impurities-Fe and La in the 
cation sublattice, and vacancies-a single missing B. For 
the first time, through the use of explicitly correlated 
many-body wavefunctions in variational Monte Carlo 
(VMC) and diffusion Monte Carlo (DMC) calculations, 
we conclusively show that pure CaBg is a semiconductor 
with an Y-point gap of - 1.3 eV. Using HF and DFT 
methods, we have shown that i) La is not essential for 
ferromagnetism in the system, while ii) Fe substitutional 
impurities give rise to ferromagnetism, iii) B vacancy, 
and Fe and La substitutional impurities make the system 
metallic, iv) CaBg with La substitutional impurities has 
a small electron-like Fermi surface (FS), explaining some 
of the observations in ARPES experiments 0,0]. Since 
experimental information about atomic structure around 
these defects is essentially absent, our studies should be 
understood as an exploration of various possible scenar- 
ios. Nevertheless, we believe our results suggest expla- 
nations to most of the experimental observations at least 
qualitatively, and give an accurate quantitative estimate 
of the gap in the pure system, thus clarifying the overall 
physics of CaB 6 . 

We performed our HF and DFT calculations using 
CRYSTAL98 and CRYSTAL03 codes [UH^. We used 
both all electron basis sets and ECPs in our calculations 
and they gave similar results. For QMC calculations we 
used the same basis sets with ECPs as used in HF. For 
QMC calculations, we used a Slater- Jastrow wavefunc- 
tion of the form 
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The Jastrow correlation part depends on electron-ion dis- 
tances ru and electron-electron distances and the 
term includes electron-electron cusp terms and altogether 
contains 10 variational parameters. In the first step, 
these variational parameters are optimized by minimiz- 
ing the variance in the local energy, El = H being 
the Hamiltonian of the system in which the core electrons 
are replaced by pseudopotentials (for Ca the Ne-core elec- 
trons were eliminated). Finally, a large part of the vari- 
ational bias is eliminated by fixed-node DMC. Details of 
these methods have been discussed in the literature . 

Pure CaBe: As is typical with these methods, HF gave 
a large X-point gap of ~ 5.50 eV, while the use of hybrid 
functional B3LYP reduced the gap to - 0.86 eV. PW91 
generalized gradient approximation (GGA) gave a con- 
ducting state and a small band overlap at the X-point, as 
found in all other such calculations in the literature (see 
fig.Q). For QMC calculations we used a 2 x 2 x 2 super- 
cell so that the X-point of the primitive BZ is folded on 
to the T-point of the BZ corresponding to the simulation 
cell. Since we are interested in the gap at the X-point of 
the primitive BZ, we can work with wavefunctions at the 
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FIG. 1: Band structure of pure CaB6 calculated with HF, 
B3LYP functional and PW91 GGA functional. 



the Slater part of the many-body wavefunction for our 
QMC calculations were taken from HF calculations. For 
the ground state wavefunction, the N lowest energy or- 
bitals were included. For the excited state, we promoted 
one electron from the top of the valence band to the bot- 
tom of the conduction band at the T-point of simulation 
BZ which corresponds to the folded X-point. Variational 
principle guarantees that the calculated energy is a strict 
upper bound for the lowest energy state of the symmetry 
imposed by the trial function. Since we are interested in 
direct gap, the excited state of interest in not the ground 
state of the particular symmetry (which is determined by 
the fc-point occupation in a solid). It is well known that 
in such cases the DMC energy is not necessarily a strict 
upper bound of the excited state energy. However, fixed- 
node DMC gives a good estimate of the excited state en- 
ergy. In fact, it is exact if the nodal structure of the trial 
function W(i?) is exact, and in most practical cases the 
fixed-node constraint is dominant and leads to energies 
which are higher than the exact energies for a given state 
(see ref. [24| , for example) . We optimized the variational 
parameters in the Jastrow part only for the ground state 
wavefunction, and used the same set for the excited state. 
The precise value of the gap we obtained in our DMC 
calculations is 1.3(2) eV. This is in very good agreement 
with all the experiments that give an estimate of the bulk 
band-gap [13, [ljj [lj| . Obviously, this result is based on 
fixed-node DMC single determinant trial wavefunctions. 

Our results based on the fixed-node single-determinant 
trial wave function DMC method, which is a many-body 
approach different from previous calculations, suggest 
that the pure CaBg compound is a narrow-gap semicon- 
ductor. Hence we conjecture that the metallic character 
seen in transport measurements have originated from de- 
fects in the system. This result also suggests that the 
theories with unconventional mechanisms explaining FM 
in this system using early DFT band structures might 
not be correct and we believe that our calculations sup- 
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TABLE I: Defect formation energies in eV calculated with HF 
and B3LYP. 



Defect 


HF 


B3LYP 


Fe-substitution 


15.6 


6.9 


B- vacancy 


15.8 


14.8 


La-substitution 


-0.5 


-2.8 



follow with calculations of CaBg with defects, proposing 
explanations for the origin of FM and metallic behavior 
as observed. As we mentioned, experimental information 
about atomic structures of possible defects in hexaborides 
is not available. Therefore our calculations of the hex- 
aboride with defects only serve as models to illustrate 
the qualitative impact of imperfections on the electronic 
structure. 

CaB e with defects: We present here results from our 
calculations of three different defects in CaBg. We con- 
sider substitution of Ca in the lattice by Fe and La, and 
an isolated B vacancy. The first one is motivated by the 
suggestions that magnetic ions, possibly Fe or Ni, are re- 
sponsible for FM in this system. Although there have 
been indications that Fe and Ni reside close to the sur- 
face layer [2(| , there are also experiments suggesting that 
they penetrate substantially into the bulk |25j • We treat 
these as bulk substitutional impurities to avoid compli- 
cations of surface effects such as relaxation and recon- 
struction in our calculations. It has also been speculated 
that FM with high T c is due to boride formation-FeB or 
Fe2B-which have T c 's of 598 K and 1015 K respectively. 
However, this has been contradicted by other works |18j . 
Since the experimental evidence is not conclusive, we do 
not consider any possible boride formation either. B va- 
cancy would presumably be the most natural defect in 
a nominally B deficient sample. Although no FM was 
seen in B deficient samples in ref. > we still study this 
system to understand its electronic properties. La sub- 
stitution is motivated by the original observation of FM 
in the La substituted compound Cai_ x La x B6. 

Concentrations of various defects in conditions of ther- 
mal equilibrium depend on their formation energies (FE). 
The defect formation energy is defined in the usual way, 
as the difference between the binding energies of the pure 
solid and the solid with defect. Since there are no hints 
about structural relaxations around defects from experi- 
ments we did not consider such relaxations. At least in 
some cases these have been shown to be negligible [l7j . 
However, such relaxations cannot be ruled out a priori 
and therefore our formation energies (Tabled) are most 
probably biased towards the defect formation being less 
favorable. Nevertheless, since we are interested in a ba- 
sic, qualitative understanding of the electronic structure 
changes in presence of defects, we believe that our model 
calculations are sufficient for this purpose. 
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FIG. 2: (Color online) Band structure of CaB6 with 
3.7% Fe impurity with HF, B3LYP and PW91 functionals. 
Black(solid) and red(dotted) lines represent up and down spin 
energies in the FM ground state. Fermi energy in all cases has 
been set to 0. 

from impurities in boron, or by some other means during 
sample preparation, would be present in minute amounts 
because of its relatively large FE. The smallest impurity 
concentration we have been able to study is ~ 3.7% of Fe 
doping (one Ca replaced in a 3 x 3 x 3 supercell) because 
of system size limitations. Going down from 12.5 % to 
3.7 %, an Fe impurity always stabilizes a FM state over 
a paramagnetic one. As an origin of this FM state, we 
note that there is a splitting of Fe d states due to the 
ligand field of the B octahedra. In presence of this cubic 
field, the five degenerate d states split into two different 
energy levels: triply degenerate t2 S states, and doubly 
degenerate e g states. The result is a | spin on Fe with 
a — i spin over surrounding B atoms that forms a local 
magnetic moment. A FM state will form when these local 
moments interact either through hybridized Fe-B states, 
or in regions where there is a high Fe concentration, such 
as near the surface. 

Fig. shows the band structure plots for a CaB6 sys- 
tem with a Fe impurity, calculated with HF, standard 
B3LYP and PW91 functionals. From the band struc- 
ture, the system is seen to be metallic at the B3LYP level. 
From the position of the Fermi level it is also clear that 
the system would have a small hole-like FS. We consider 
these results qualitative in nature since the concentration 
of the iron defects is artificially high due to the size of 
the supercell. Nevertheless, the presence of the Fe lo- 
calized moment will be similar in smaller concentrations 
since the t^g states are quite deep and show very little hy- 
bridization with the neighboring boron states. The most 
important changes one can expect might be actually in 
the minority spin-channel where larger supercells might 
enable better state relaxations with possibly larger spin- 
density compensation resulting in overall smaller spin po- 
larization. 

Next we consider a B vacancy in a 2 x 2 x 2 system 
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FIG. 3: B3LYP band structure and density of states of 3.7% 
La-doped CaB6 . The Fermi energy is indicated by the dashed 
line. 

moving a B has the effect of taking away aBp valence 
electron from top of the valence band. The system be- 
comes metallic and the position of the Fermi energy tells 
that it has a hole-like FS at the T-point of the BZ of the 
2x2x2 cell. Although our calculations show that a 
magnetically ordered state is favored, the energy differ- 
ence with a paramagnetic state is only ~ 0.05 eV per 
primitive cell. This energy difference is too small to be 
taken as an indication of a magnetic ground state, and 
in fact, can be an artifact of the B3LYP method. This 
conjecture is supported by the GGA results in ref. 
where B vacancy is not seen to give rise to any apprecia- 
ble magnetic moment. 

Now we come to our study of La substitution of 
Ca. Because of the initial observation of magnetism in 
Cai-^La^Bg, it was thought that La plays a role in the 
origin of FM in the system. This makes the study of La 
substitution particularly interesting. Also note that the 
formation energy of a La substitutional impurity turns 
out to be negative indicating that any La atoms present 
during sample preparation would readily get incorpo- 
rated into CaBg. In our HF calculations we find that 
La-doped CaBg, down to an impurity concentration of 
1.6%, has a paramagnetic ground state. This is found to 
be the case also at the B3LYP level of approximation. In 
any case, since HF is known to favor a high-spin state, 
one can reasonably expect the La-doped system to be 
nonmagnetic. This, along with our prediction of FM in 
Fe doped CaBg, lends support to the claims that FM in 
this system is an extrinsic effect due to some magnetic 
atoms, and is not related to the presence of La. In Fig.|3 
we show the band structure and density of states for a 
system with 3.7% La-doping. The system is found to 
be metallic. From the position of the Fermi surface we 
find that the system has a small electron-like FS at the 
r-point in the BZ. Comparing band structure in Fig. 
with that of a (3 x 3 x 3) supercell (not shown) of pure 

/~< t-» p • _ 1.1- n . r i! . ti. 



energy with only minor changes in energy bands. Thus 
La replacing a Ca atom acts as an electron donor in the 
system. We believe that the electron-like Fermi surface 
seen in ARPES measurements 00 might have origi- 
nated from La substitutional impurities. 

Finally we study the effects of a B vacancy and a Fe 
substitutional impurity being present together in CaBg. 
This is motivated by the observation in ref. 0] that al- 
though CaBg, CaBg + 5 and Cai-^La^Bg prepared with 
lower purity B exhibit FM, no such effect is seen in 
CaBg_5. This suggests that although Fe substitutional 
impurities give rise to FM, B vacancies neutralize that 
effect and when both are present in the system it is no 
longer FM. This is precisely borne out by our calcula- 
tions. We have studied two extreme scenarios as far 
as the positions of the B vacancy and Fe impurity are 
concerned-i) they are on neighboring sites, and ii) they 
are as far apart as possible in our 3x3x3 simulation 
cell. In both cases a nonmagnetic state is preferred over 
a FM state in our HF calculations. Again, HF typically 
preferring a high-spin state, we can take this as a reason- 
able indication that a B vacancy destroys the FM state 
created by a Fe impurity. 

In conclusion, using explicitly correlated many-body 
wavefunctions, we have conclusively shown that pure 
CaBg is a narrow-gap semiconductor with an A-point 
gap of 1.3(2) eV. Our HF and DFT calculations show 
that possible imperfections in the crystal-B vacancy, Fe 
or La substitutional impurities-drive the system metallic, 
explaining the apparent anomalies in the experimental 
results. Out of these, La substitutional impurities give 
rise to a small electron-like Fermi surface, which was seen 
in experiments but had no clear explanation. We have 
also shown that FM in this system is related to magnetic 
impurities, e.g., Fe and probably has no relation with the 
presence of La. Another interesting fact emerging from 
our calculations, explaining the fact that FM was never 
seen in CaBg_,5 samples in Ref. [l^, is that a B vacancy 
destroys any magnetic moment created by Fe impurities. 
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